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qq ' Abstract. A probability density characterization of multipartite entanglement is 

CN ■ tested on the one-dimensional quantum Ising model in a transverse field. The average 

and second moment of the probability distribution are numerically shown to be good 
indicators of the quantum phase transition. We comment on multipartite entanglement 
00 i generation at a quantum phase transition. 
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Quantum phase transitions are characterized by nonanalyticity in the properties of the 
states of a physical system [1]. They differ from classical phase transitions in that they 
occur at zero temperature and are therefore driven by quantum (rather than thermal) 
^ | fluctuations. 

The research of the last few years has unearthed remarkable links between quantum 
phase transitions (QPTs) and entanglement [2, 3, 4, 5]. The study of these inherently 
quantum phenomena has mainly focused on bipartite entanglement, by using the entropy 
of entanglement [6], i.e. the von Neumann entropy of one part of the total system in 
the ground state. Notwithstanding the large amount of knowledge accumulated, the 
properties of the multipartite entanglement of the ground state at the critical points 
of a QPT are not clear yet. This is also due to the lack of a unique definition of 
multipartite entanglement [7]. Different definitions tend indeed to focus on different 
aspects of the problem, capturing different features of the phenomenon [8] , that do not 
necessarily agree with each other. This is basically due to the fact that, as the size of 
the system increases, the number of measures (i.e. real numbers) needed to quantify 
multipartite entanglement grows exponentially. For all these reasons, the quantification 
of multipartite entanglement is an open and very challenging problem. 

In the study of a QPT the above-mentioned problems are of great importance. The 
evaluation of entanglement measures bears serious computational difficulties, because 
the ground states involve exponentially many coefficients. The issue is therefore to 
understand how to characterize entanglement, e.g. by identifying one key property that 
can summarize its multipartite features. Our strategy will be to look at the probability 
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density function of the purity of a subsystem over all bipartitions of the total system. 
The average of this function will determine the amount of global entanglement in the 
system, while the variance will measure how well such entanglement is distributed: a 
smaller variance will correspond to a larger insensitivity to the choice of the bipartition 
and, therefore, will witness if entanglement is really multipartite. 

This approach, introduced in [9] , makes use of statistical information on the state 
and extends in a natural way the techniques used for the bipartite entanglement. It is 
interesting to notice that the idea that complicated phenomena cannot be "summarized" 
by a single (or a few) number(s) was already proposed in the context of complex systems 
[10] and has been also considered in relation to quantum entanglement [11]. We applied 
our characterization of multipartite entanglement to a large class of random states 
[12, 13], obtaining sensible results [9, 14]. 

In this article we will characterize in a similar way the multipartite entanglement 
of the (finite) Ising model in a transverse field. Our numerical results will corroborate 
previous findings and yield new details about the structure of quantum correlations near 
the quantum critical point. 

2. Probability density function characterization of multipartite 
entanglement 

We shall focus on a collection of n qubits and consider a partition in two subsystems 
A and B, made up of n& and % qubits (jia + tib = n), respectively. For definiteness 
we assume ua < fis- The total Hilbert space is the tensor product TC = Ha <8> T~Cb with 
dimensions dimH A = N A = 2 n - 4 , dimH B = N B = 2 nB and dim H = N = N A N B = 2™. 
We shall consider pure states 

N-l N A -1N B -1 

m = X>i*> = E E ^fetfA> ® ik>. (!) 

k=0 j A =0 l B =0 

where the last expression is adapted to the bipartition: \k) = \ja) <8> \Ib), with a bijection 
between k and (j'a,^b)- Think for example of the binary expression of an integer k in 
terms of the binary expression of (ja,Ib)- We define the purity (linear entropy) of the 
subsystem 

kab(\^)) = Tr A p^, p A = Tr B p, p=\^){ip\, (2) 

Tta (Tr#) being the partial trace over subsystem A (B), and take as a measure of the 
bipartite entanglement between A and B the participation number 

Nab = vr^, (3) 

that measures the effective rank of the matrix pa, namely the effective Schmidt number 
[15]. The quantity uab = log 2 Nab represents the effective number of entangled qubits, 
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given the bipartition (pictorially, the number of bipartite entanglement "links" that are 
"severed" when the system is bipartitioned). By plugging (1) into (2) one gets 

/N A -1 N B -1 \ ~ X 

Nab(\iI>)) = I E z jiZj'iZj'i>Zji> I • (4) 

\j,j'=0 l,l'=0 J 

This is the key formula of our numerical investigation. 

Clearly, the quantity Nab will depend on the bipartition, as in general entanglement 
will be distributed in a different way among all possible bipartitions. We are pursuing 
the idea that the density function p(Nab) of Nab yields information about multipartite 
entanglement [9]. We note that 

\<N AB = N ba < N A (< N B ), (5) 

where the maximum (minimum) value is obtained for a completely mixed (pure) state 
Pa- Therefore, a larger value of Nab corresponds to a more entangled bipartition 
(A,B). Incidentally, we notice that the maximum possible bipartite entanglement 
Nab* = N™ 3 * = 2'"/ 2 ' can be attained only for a balanced bipartition, i.e. when 
ha = [n/2] (and ns = [(n + l)/2]), where [x] is the integer part of the real x, that 
is the largest integer not exceeding x. We emphasize that the use of the inverse purity 
(linear entropy) (3) is only motivated by simplicity. Any other measure of bipartite 
entanglement, such as the entropy (or any Tsallis entropy [16]) would yield similar 
results. 



3. Entanglement distribution: critical Ising chain in a transverse field. 

We now apply the characterization of multipartite entanglement to the quantum Ising 
chain in a transverse field, described by the Hamiltonian 

n—l n n 

H = -9J2 «+i - (! - 9) E of + e < (6) 

i=l i=l i=l 

(with open boundary conditions, a being the Pauli matrices). Notice that we added a 
(small, site independent) longitudinal field e. If e = 0, it is known from conformal field 
theory [17] and numerical simulations based on accurate analytical expressions [3] that 
at the critical point g = g c = 1/2 the entanglement entropy 

Sab = - Tr A (pyi log 2 Pa) (7) 
diverges with a logarithmic law 

S A B~-xlog 2 £. (8) 
o 

Here entanglement is evaluated by considering a block A of contiguous spins whose 
length I is less than one half the total length n of the chain. Due to (approximate) 
translation invariance, in our approach this is equivalent to considering the average 
entanglement over a subset of the bipartitions of the system (that tend to be balanced 
when i tends to n/2). 
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3.1. A typical distribution 

We intend to evaluate the distribution of bipartite entanglement over all balanced 
bipartitions and, therefore, the multipartite entanglement. Here and in the whole article, 
the Hamiltonian will be exactly diagonalized in order to obtain the ground state, then 
Nab will be explicitly evaluated as a function of g and its distribution plotted. The 
results are exact, but the quantum simulation time consuming and for this reason n 
cannot be too large. 

The distribution of the participation number Nab as g varies, for n = 10 qubits and 
e = 0, is shown in Fig. 1. We notice that the distribution is always well-behaved and 
bell-shaped, being practically a 5 function for g < 0.1 and g > 0.75. For this reason, 
one can get a satisfactory characterization of multipartite entanglement by looking at 
its mean value and width 

li= {Nab), a 2 = ((N A B-ti) 2 ), (9) 

where the average (• • •) is evaluated over all balanced bipartitions. We recall that /i 
defines the amount of entanglement while the inverse width a -1 describes how fairly 
such entanglement is distributed. We notice that the width a is maximum at g = 0.5, 
while the average entanglement /i is maximum at g = 0.56. Observe that no singularities 
can be expected for a number of spins as small as n = 10, yet the behavior of both 
quantities clearly foreruns the quantum phase transition at g = g c = 1/2. In this sense, 
both a and [i appear to be good indicators of the QPT. 

3.2. Average and width 

Let us consider the full Hamiltonian (6) when the longitudinal perturbing field is small. 
In Fig. 2 we plot fi and a, respectively, vs g for the ground state of the Hamiltonian 
(6), when n = 9, for different values of e (ranging from to 10~ 2 ). We notice a 
very different behavior of the two quantities. The average /i is very sensitive to the 
longitudinal perturbation. In the region g ~ 1, where the ground state is approximately 
a GHZ state when e = 0, the average entanglement is strongly reduced even for a very 
small value of e (~ 10 -6 ). This is basically due to the fact that the superposition 
| all spins up) + | all spins down) (yielding \i = 2) is very fragile and the ground state 
collapses in one of the two (degenerate) classical ground states (yielding \i = 1), the Z 2 
symmetry being broken. On the other hand, near the maximum, \x is more robust and a 
larger perturbation (e = 10~ 2 ) is required to counter larger values of (1 — g) and modify 
the behavior of fi. 

The behavior of a is different. When e < 10 -2 the curves are not modified by the 
presence of the longitudinal field. This is due to the fact that in the region where \x is 
reduced by the presence of e, a is already near to (a GHZ state has a = because it 
is invariant for permutation of the qubits, see [9]). Of course, a sufficiently large value 
of e affects also a, reducing it (but not modifying the shape of o~(g)). 
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Figure 1. Distribution function of the participation number Nab over all balanced 
bipartitions for the Hamiltonian (6) when e = and n = 10. The distribution is 
always bell-shaped. Its width is maximum at g = 0.5, while its average entanglement 
(indicated by a black arrow) is maximum at g = 0.56. Notice the different scales on 

the ordinates. The number of balanced bipartitions is n v = I , ) = | ) = 252. 

\[n/2\ 
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Figure 2. (a) Average fi and (b) standard deviation a of Nab over all balanced 
bipartitions for n — 9 sites in the 1-D quantum Ising model in a transverse field of 
strength 1 — g and a small longitudinal field of strength e. Squares: e = 0; stars: 
e = 1CP 6 ; triangles: e = 10~ 4 ; diamonds: e — 10~ 2 . 



We shall now focus on the critical region. It would be tempting to take a small 
value of e (say, e = 10~ 6 ) in order to get rid of the spurious residual entanglement at 
g — 1 (and obtain a bell-shaped function for fi — as well as for a). However, since we 
aim at a precise determination of the coordinates of the maximum, which is unaffected 
by small values of e, we decided to work with e = 0. 

3.3. Purely transverse Ising chain 

In Fig. 3 we evaluate the average and standard deviation for e = (purely transverse) 
Ising chains of increasing size (from 7 to 11 sites). In Fig. 3(a) we distinguish different 
zones. For g = the ground state (gs) is factorized and \i — 1. If g ~ 1 the gs is 
approximately a GHZ state (a combination of the gs's of the classical Hamiltonian) . 
The most interesting region is around the value g = 0.5, where for an increasing number 
of qubits there is a more and more pronounced peak of \x. This is in qualitative agreement 
with other results obtained using the entropy of entanglement. 

The width of the distribution of Nab versus g is shown in Fig. 3(b). We will 
comment later on the behavior of this quantity, that yields useful additional information 
about the structure and generation of multipartite entanglement (information that 
would not be available for an entanglement measure constituted by a single number). 
Also in this case we can distinguish several regions in the plot. Moreover, the coupling 
g corresponding to the peak of a (that we denote <r max ), does not coincide with that 
corresponding to the peak of \i (that we denote /i max ) : 

g (Omax) < g (/imax) • (10) 

In other words, for a finite spin chain, the width of the distribution is not maximum 
when the amount of entanglement is maximum. 

We notice that, by increasing n, both maxima are shifted towards the center of the 
plot g — ► g c = 0.5. In Fig. 4(a) we plot the values of the coupling constant g at fi max 
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Figure 3. (a) Average fi and (b) standard deviation a of Nab over all balanced 
bipartitions (from n = 7 to 11 sites) for the purely transverse 1-D ising chain. Full 
squares: 11 sites; open triangles: 10 sites; open squares: 9 sites; full triangles: 8 sites; 
open diamonds: 7 sites, [i can be viewed as a measure of the average multipartite 
entanglement, while cr -1 can be viewed as a measure of how fairly this entanglement 
is shared. Both fi and a are good indicators of the QPT that takes place at g = 0.5. 
Interestingly, cr max precedes /z max . 



versus the number of sites n. The numerical result can be fitted with the (arbitrary) 
function 

<K/w) = 0.5 + 2 ^ qn 5 Q 43 ™ °' 5 = 9c (11) 

rr + 3.09 n — 35.59 

The plot of g(cr max ) versus n is shown in Fig. 4(b), the fit being 

Q ((Tmax) = 0.5 + — — ™ 0.5 = g c . (12) 

yv ; n 2 - 13.01 n + 46.39 y V ; 

Notice that the fit (11) is very accurate, while (12) is valid within one standard deviation 

(namely a few percent), as can be seen in Fig. 4. From Fig. 4 and Eqs. (11)-(12) one can 

argue that the amount of entanglement (the mean of the distribution) and the maximum 

width of the distribution of bipartite entanglement can detect, in the limit of large n, 

the QPT. 

We shall henceforth focus on /i max and cr(/i max ) = cr(g(/x max )) (the value of cr when 
the amount of entanglement is maximum), rather than <7 max (whose behavior is anyway 
similar). In Fig. 5 we plot these quantities vs the number of spins n. They are fitted by 
(for n > 6) 

/imax =2 + 0.019 (n- 6) + 0.007 (n-6) 2 , (13) 



cr(/i max ) = - 0.077 + 0.11 VrT^6. (14) 
We also evaluate the relative width at maximum entanglement 

frel = Cr(At max )/A t ma X , (15) 

shown in Fig. 6, that will be useful in the following discussion. The fitting curve in Fig. 
6 is not independent, but is rather derived from Eqs. (13)-(14): 



-0.077 + 0.11 
2 + 0.019 {n - 6) + 0.007 (n - 6)' 



°"rcl = ~ ; „ , ^ . —, ^2- (16) 
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Figure 4. Coupling constant g corresponding to (a) // max and (b) <7 max versus n. 
Notice that g(cr max ) < <?(^ max ) (at fixed n) and that already for small n(= 7), g(<T m ax) 
differs from g c = 1/2 only by a few percent. The error bars (one standard deviation) 
are explicitly indicated. 
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Figure 5. (a) Entanglement /Lt max and (b) standard deviation at the maximum 
entanglement cr(/^ m ax) vs n. The error bars (one standard deviation) are explicitly 
indicated. 



4. Discussion 

Both fits (13)-(14) imply that the entanglement indicators a and fi diverge with n at 
the QPT. This conclusion is particularly significant: the amount of entanglement goes 
to infinity but so does the width of the entanglement distribution. In particular, this 
leads to two possible scenarios, depending on the behavior of a rc \ defined in (15): 

(i) er rc i n ^2$ q T n case divergence of /i max is stronger than that of cr(/x max ) ~ 
cr max . This means that at the QPT the entanglement of the ground state is 
macroscopically insensitive to the choice of the bipartition. Accordingly, the QPT 
yields a fair distribution of bipartite entanglement and is therefore a good tool for 
generating multipartite entanglement. This conclusion could pave the way towards 
a deeper understanding of the relation among entanglement, QPTs and chaotic 
systems (that are known to generate large amounts of entanglement [13, 18]). 




(ii) <r rel — > c>0 (eventually oo). This situation would have profound consequences on 
our comprehension of the relation between a QPT and the generation of multipartite 
entanglement. In particular, the strong divergence of cx(/i max ) (of order equal to or 
larger than that of /i max ) would imply that the distribution of entanglement is 
not optimal, inasmuch as it is not fairly shared. This means that the amount of 
entanglement of non-contiguous spins partitions macroscopically differs from that 
of contiguous ones. 

Our results, although not conclusive due to the relatively small value of n reached in our 
numerical analysis, appear to indicate that (i) is the most probable scenario: indeed, 
from Eq. (16), that in turn is a consequence of Eqs. (13)-(14), we infer that for large n 

o-rd ~ rT 3 / 2 . (17) 

In general, if one assumes that the behavior of /i max and cr(/x max ) vs n (and in particular 
the convexity of the two curves) does not change for larger n, one can conclude that a rc i 
vanishes for n — > oo. 

Another important observation, related to the "entangling power" of evolutions [19], 
is the following. Although our numerical results seem to favor the first scenario, namely 
a well distributed multipartite entanglement generated by the quantum phase transition, 
such entanglement is not so large. Indeed, a typical n-qubit state is characterized by [9] 

fi oc 2 n/2 , a = const, (18) 

namely an exponentially large amount of entanglement, that is also very well distributed. 
These typical states are efficiently produced by a chaotic dynamics [13, 18]. In general, 
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one observes a very rapid growth of the (effective) Schmidt number (3) at the onset 
of chaos and for all these reasons, quantum chaos is a much better multipartite 
entanglement generator than a critical Ising chain. This conclusion seems to be valid for 
other spin Hamiltonians as well. Notice that the entangling power (and/or entanglement 
generation) of a QPT is better compared to that of a chaotic system [18] (in that they 
are both obtained by varying one or more coupling constants), rather than that of a 
quantum evolution [19]. On the other hand, unlike in a chaotic system, in a QPT one 
focuses on the features of the ground state. 

The entanglement generation at a QPT and the physical features of this 
entanglement [20, 21] deserve additional investigations. The participation number or 
the entropy of entanglement (or any other sensible measure) are related to the global 
structure of the state. It is therefore reasonable to expect that many observables might 
be necessary in order to characterize multipartite entanglement. The approach we 
propose [9, 14], based on the calculation of the probability density function of bipartite 
entanglement, has the advantage of making use of statistical information on the state 
of the system and characterizes multipartite entanglement by extending techniques that 
are widely used in the analysis of its bipartite aspects. We have seen that when the 
density functions are well behaved and bell-shaped, the average and second moment of 
the distribution are good indicators of the quantum phase transition. These conclusions 
must be corroborated by the study of other systems and models displaying quantum 
phase transitions, as well as by the analysis of more complex systems [10, 11]. Work is 
in progress in this direction. 
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